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Abstract 

Objectives: Mathematical models can be helpful to understand the complex 
dynamics of human immunodeficiency virus infection within a host. Most of work 
has studied the interactions of host responses and virus in the presence of active 
cytotoxic immune cells, which decay to zero when there is no virus. However, 
recent research highlights that cytotoxic immune cells can be inactive but never 
be depleted. 

Methods: We propose a mathematical model to investigate the human immunode- 
ficiency virus dynamics in the presence of both active and inactive cytotoxic immune 
cells within a host. We explore the impact of the immune responses on the dynamics 
of human immunodeficiency virus infection under different disease stages. 
Results: Standard mathematical and numerical analyses are presented for this new 
model. Specifically, the basic reproduction number is computed and local and 
global stability analyses are discussed. 

Conclusion: Our results can give helpful insights when designing more effective 
drug schedules in the presence of active and inactive immune responses. 



1. Introduction 

Human Immunodeficiency Virus (HIV) is the path- 
ogen which is responsible for the Acquired Immuno- 
deficiency Syndrome, AIDS. Once HIV enters into the 
body is detected by macrophages and a range of antigen 
presenting cells; in particular, dendritic cells and 



macrophages search, phagocytize and analyze the anti- 
gen to be presented to inactivate CD4 T-cells in a pro- 
cess called immune synapse or simply antigen 
presentation [1]. Dendritic cells have certain lectins in 
their membranes to which viral particles bind with large 
affinity, this traps HIV and during the antigen presen- 
tation process CD4 T-cells are more likely to get 
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infected by the surrounding viral particles attached to 
the dendritic cell membrane [1,2]. 

HIV mostly targets the immune response system, 
particularly, CD4 T-cells, where the viral RNA is con- 
verted into viral DNA so that when these cells are 
activated, they can replicate virus. CD4 T-cells are 
responsible for signaling other immune responses such 
as the cytotoxic and humoral responses that are invaders 
that have to be fought. In particular, Cytotoxic T 
Lymphocytes and CD8 cells direct cytotoxic immune 
response. This cytotoxic response act by lysing infected 
cells causing them to explode, thus cytotoxic cells 
remove infected CD4 T-cells from the body at a con- 
stant rate, but these cells do not directly target free virus. 
Over the time, HIV is able to deplete the population of 
CD4 T-cells in such way that the Cytotoxic response is 
never depleted. 

There have been remarkable attentions on the HIV 
dynamics among scientists in medicine and mathemat- 
ical biology. Significant efforts have been made in order 
to understand and to characterize the underlying mech- 
anism of the disease. Earlier mathematical framework 
has been considered to model HIV/AIDS dynamics 
focusing on the viral and CD4 T-cells dynamics [3—1 1] 
and the references therein. Recently, there has been 
some work to explore the impact of the immune 
response in the HIV dynamics [12—15]. Specifically, 
mathematical models incorporated explicitly the effect 
of the immune response and their results have been 
analyzed and investigated [12,14,15]. Most of work has 
studied the interactions of host responses and virus in 
the presence of active cytotoxic immune cells, which 
decay to zero when there is no virus. However, recent 
research highlights that cytotoxic immune cells can be 
inactive but never be depleted. 

In this manuscript, we propose a mathematical 
model, which further captures active and inactive im- 
mune cells in the previous models. Our main focus is to 
identify the role of both active and inactive immune 
responses in the HIV dynamics within a host. We 
analyze a modification of the existing model [12] that 
splits the cytotoxic immune cells compartment, M, into 
two compartments, inactive and active immune response 
cells. On one hand, this process of activation leads us to 
take into account that cytotoxic cells are always present 
in the body. On the other hand, inactive cytotoxic im- 
mune cells get activated through specific biochemical 
processes related to the presence of HIV. 

2. Materials and methods 

2.1. Previous models 

Let V(t) denote the average viral particle concentra- 
tion at time t assuming that when the initial time t = 0, 
an initial viral load Vq > 0 enters into the body. In ideal 
conditions, this initial viral load is going to be 



eliminated from the body at a constant rate c and this 
depends only on the virus ability to infect the immune 
cells. Therefore, as long as the infection has not been 
established, the viral load can be described by 

V=~cV (1) 

where V represents the rate of change in the viral con- 
centration per unit time. Ata very early stage of infec- 
tion, the virus doesn't find the proper conditions for 
successful replication, it decays exponentially as given 
in the expression, V(t) = V§e~ ct . Perelson and Nelson 
studied a similar situation to the one we just described 
but there is an unknown quantity (and to be determined) 
describing the creation of new viral particles [16]. To 
define this quantity, first, we consider that once the virus 
enters into the body and it will infect the CD4 T-cells, 
which are immunologically activated T-cells. Let 
T = T(t) be the average concentration of healthy CD4 
T-cells at time t at the constant recruitment rate a and 
death rate /j.. In this way T = <j — \lT describes the 
evolution of T without infection. In the absence of virus, 
the CD4 T-cells reaches the equilibrium level of a/fi 
cells per mm 3 . If we denote /? as the probability of a 
CD4 T-cell for to be infected by the HIV, then from the 
mass action principle, fiTV represents the average 
number of CD4 T-cells per unit of time that getting 
infected at time t. Thus the equation becomes 

t = a- (3TV - fiT. (2) 

Now, the average concentration of infected CD4 
T-cells is represented by the variable T* = T* (t), these 
are activated by the HIV aiming to replicate more viral 
particles. One remarkable aspect of this disease is that 
the destruction of infected CD4 T-cells was initially 
thought as a consequence of cytopathic action but later 
it was found other indirect destruction mechanisms 
such as induction of apoptosis through soluble viral 
proteins, secondary cellular death due to immunolog- 
ical hyper-activation, syncytia formation and progres- 
sive damage of the primary and secondary lymphoid 
organs [1,17]. Let 8 be the infected CD4 T-cells death 
rate, so 8T* is the average concentration of infected 
CD4 T-cells that die at a time t. Hence, the equation 
describing the variation of infected CD4 T cells take 
the form, 

T'=(3TV~5T' (3) 

Furthermore, let us consider that each CD4-T infec- 
ted cell produces rj new viral particles, then the equation 
(1) becomes, 

V = r]8T' -cV. (4) 

Where rj5T* represents the average number of new viral 
particles produced at the time t. So far, from these 
equations (2—4), we have constructed the following 
system of ordinary differential equations (ODEs) 
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T = <T-PTV-fiT 
f =$TV -dT' 

v= T]dr -cv. 



(5) 



The system (5) is the simpler model for HIV and 
immune system dynamics that can be considered. It has 
been widely studied by [13,18] among many other au- 
thors. Different models, which are derived from this 
system, have been mathematically analyzed in the 
literature [2,19—21]. However, this ODE model does 
not take into account the immune response (cytotoxic 
immune response or humoral immune response) effect, 
which is a key phenomenon in the dynamics of HIV 
within a host. Nowak and Bangham have included this 
effect and given as the following ODE system [12]: 



T = <r-l3TV-fiT 

f" =PTV-dT' -yT'M 

M = aT'M -5M 

V = 7]5T' -cV 



(6) 



This system includes M = M(t), which corresponds 
to the average concentration of active cytotoxic immune 
response cells, i.e. cells capable of eliminating infected 
CD4 T cells through cytotoxic action with a probability 
y. Then, the term yT'M denotes the average concen- 
tration of infected CD4 T-cells that are destroyed. 



number of inactive immune response cells that become 
active. The natural death rate of both inactive and active 
immune cells is denoted by p. Hence, the mathematical 
model can be written as: 

dT 

-=a-(3TV-,T 



dT' 
dt 

dM _ 
dt 

dM' 
dt 

dV_ 
dt 



= I3TV~ST' -yT'M 
■X- \f/T"M - pM 
= aT'M' +\j/T'M- pM' 
r)T* - cV 



(7) 



Note that last equation for system (6) differs from the 
corresponding equation on system (7)in the term rjT*. 
Within the system (7), the factor 5 is excluded because 
we want to consider all the cellular production of virus, 
not only those who are released by the cell when it dies 
out. In this model we neglect, as in [7,12,16], the loss of 
virus during the infection. Once the model is formulated, 
standard mathematical analysis is carried out including 
the stability analysis based on the basic reproduction 
number R 0 (shown in Appendix). 



2.2. HIV model with both active and inactive 
immune cells 

We propose a new model for HIV infection in the 
presence of both active and inactive immune cells, 
which correspond to cytotoxic immune response of the 
body by including a class of non-active immune 
response cells to avoid the extinction of immune 
response in the absence of HIV and so it never dies out. 
The population is divided as follows: T = T{t) and 
T* = T*(f) describe the average concentration of 
healthy and infected CD4 T-cells, respectively. 
M = M(t) and M* = M*(t) correspond to the average 
concentration of inactive and active cytotoxic immune 
response cells. V = V(t) is the average concentration of 
viral particles at time t. The time variation for T, T* and 
V are modeled in similar manners as the ones described 
in the introduction. The active cytotoxic immune 
response cells, M*, kill infected cells by an average 
quantity yT*M*. The detection of infection in the body 
produces aT*M* active cytotoxic immune cells, there- 
fore, a is a cytotoxic immune response activation rate. It 
is biologically meaningful to consider y > a because it 
implies that immune response cells kill more cells than 
they replicate themselves by this process. The inactive 
immune response is self-produced at a constant rate of A. 
The infected cells stimulate the inactive immune 
response cells at a rate of \p. Hence \pT*M indicate the 



3. Results 

In this section, a selection of the numerical simula- 
tions from the model (7) is presented. The parameter 
values that were used to obtain these simulations are 
given in Table 1. Two situations have been considered 
according to the changes in parameters (3, namely, 
R 0 > 1 and R 0 < 1. First, in the case of R 0 > 1, Figure 1 
illustrates the HIV dynamics of a patient at an initial 
stage of infection with no drug treatment under four 
distinct values of (3. It is assumed the patient has a 
normal count of CD4 + T-cells around 1000 cells per 
unit volume at the beginning. As we can see in the last 
panel of Figure 1 , it is typical to observe the viral par- 
ticles reaches the peak at around between day 20 and 40, 
then, they decrease and remain low for the rest of the 
simulation time duration. Also, it is clear that as greater 
the value of /3, greater the average of infected CD4 T 
cells and viral particles. Also, the peak timing of 
infected CD4 T cells and viral particles occurs earlier as 
the value of /? gets larger. 

Once the entrance of virus is perceived, the organism 
responds by activating immune response cells that fight 
against the virus as shown in the two middle panels of 
Figure 1. As the dynamics of the infection is explored 
using four distinct values of 13, each one of results show 
different levels of intensity for active and inactive 
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Table 1. 


TlPTiniti r»n r\T ncir'Qrnptpi'c *anrl iraliiAG 
J_^C11111L1U11 Ul VJaL til 1 1L IC1 S dllU V dlUCs 






Param. 


Description 


Value 


Reference 


Tr. 

1 0 


Tn iti al vn Iiip Ynv nni nfpptpn f^Xl^— 1— i pp 1 1 c 

111111(11 V £11 Ut l\Jl Ullllllt,l_- ICU V^L/t 1 i 


1000 




T* 

-'o 


Initial value for infected CD44-T cells 


0 




1V1Q 


Initial valnp for N^n-apti vp ittittiiitip ppIIq 

lllllltll VulUt l\Jl I>Ull Uvll 1 1 1 11 1 J LI J L 1 1 1 .i 


0 




Ml 


Initial value for active immune cells 


1 




v„ 

'0 


Initial walnp Tr\v vinic 
1111 I Id 1 V dlUC 1U1 V 11 LIB 


0 01 




/j 


Slnnrpp tpnn for nmnfprtpH f^D4 X ppIIq 

LJ \_J LI A V w I ^ 1 111 L\Jl Lll 11111 L< w LVu V IV i ± ^.^113 


10 mm 3 H _1 

1 \J 111111 \Jl 


T221 


P 


T? atp (^TiZl T ppll npprvmpc iTTrpr*tpri r\\? \71n1c 
JVtlLC v---L* t T 1 LC11 UCLUlllCa 1111 LCU. Uy V 11 Lla 


1 S*1 0 -5 mm 3 H _1 
Z,. J lu 111111 u 


L ZZ J 


il 


Death rate of uninfected CD4 T cell 


0.01 d _1 


T221 


1 


Niimhpr nf vinic n,mniippri nv ppI 1^ Ivciq 

INU.lllL'tl Ul VH pi UUUCtU UV VjGIIO lyolo 


500 


T221 


,5 


l~)patVi ratp of infpptpH r^T~)4 X ppIIs 


0.26 d~' 


T221 


c 


Clearance rate of virus 


2,4 d" 1 


[22] 


a 


Rate of immune response proliferation 


5*10~ 5 mmV 


[22] 


P 


Death rate of immune response 


0.1 mmV 


[22] 


7 


Rate actively infected cells deleted by CTL 


2*10~ 3 mm 3 r' 


[22] 


A 


Source term for immune response 


5 mmV 


[22] 




Rate of immune response activation 


2*10~ 3 mm 3 r' 





immune responses to help reduce infected virus cells. 
Next, Figure 2 displays the dynamics of infection of a 
patient at an initial stage of infection when R 0 < 1. 
Reduction in (3 (below a certain thresh hold number) 
leads to the reduction on the initial outbreak of infection. 
Therefore, when R 0 < 1, infection would not establish in 
the body (the number of infected virus goes to zero 
shown in the last panel of Figure 2). More extensive and 
rigorous mathematical analysis and numerical simula- 
tions will be needed for clarifying the impact of active 
and inactive immune responses on the HIV dynamics. 

4. Discussion 

A primary focus of this paper has been on mathe- 
matical modeling of the HIV infection dynamics within 
a host in the presence of both active and inactive im- 
mune responses. The present study reviewed the previ- 
ous mathematical models from the simple to the model 
with immune responses. We propose a new model, 



which is a natural extension of the existing model by 
taking account the immune response that should remain 
always in the body. 

More rigorous mathematical analysis should be and 
will be done in our future research. An initial stage of 
the HIV dynamics is taken as a numerical example. 
When one is at an advanced stage of the disease, the 
HIV dynamics would be different from what we have 
shown in this study. Therefore, various scenarios will be 
further need to be investigated. 

This work can suggest that more successful strategies 
to control HIV infection when the immune system is 
more efficient in eliminating infected cells. If the im- 
mune system is efficient (strong immune response), then 
it is possible to maintain the effectiveness of treatment 
for a longer time and at lower. Such strategies should 
contribute to build effective immune memory against 
the virus, and hence reduction of treatment costs and the 
side effects associated with the patient on multiple HIV 
therapies. The further study, which involves optimal 
drug interventions, will be carried out in our future 




Figure 1. Parameter /3 was varied to illustrate its effect on the infection's evolution at (3 = 0.000015 (dotted), /3 = 0.000025 
(line-dot); j3 = 0.000035 (dashed) and /? = 0.000045 (solid) when R 0 > 1. 
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Figure 2. Parameter /3 = 0.000004 is used to illustrate no infection occurs when R 0 < 1 . 



research. It is certainly known that HIV dynamics is 
more complicated than the one reflected by this given 
model but the findings in this paper exhibit the different 
possibilities and approaches of using mathematical 
modeling and numerical simulations in order to achieve 
new insights in more effective drug treatments. 



Proposition 1. If R 0 > 1, there exists a unique 
endemic equilibrium of the system (7). 

Proposition 2. If i? 0 < 1> the virus free equilibrium of 
the system (7) is globally asymptotically stable. 

The proofs of the above propositions and rigorous 
mathematical analysis will be given in our future research. 



Appendix. Mathematical analysis 

Lemma 1. The region Q= |(r, T* ,M,M" , 
V)eR 5 ,:T < T + T* + M + M* < ^— — V < 

[l E 

^ a — — -} is positively invariant for the system (7). 
ce i 

This result is important from the physiological point 
of view due to the fact that it ensures our mathematical 
model is biologically relevant. The positive invariance 
guarantees none of the populations, within the model, 
either goes below zero or growth without a limit. In 
other words, it assures that all populations are nonneg- 
ative and finite at every time t. 

The virus free equilibrium of the system (7) is given 
by (To.oMofi,o) = (^/ 1^,0, X/ p, 0,0) and it belongs al- 
ways to O. This equilibrium represents the situation 
without infection in the body. 

The basic reproduction number is given by the 
expression [23]: 

It denotes the number of secondary cases produced 
by an infected cell during its lifespan into a susceptible 
healthy CD4 T-cells population. Infected cells produce ? 
virions during its lifespan. These virions infect W To on 
the whole healthy CD4 T-cells population. This 
dimensionless parameter is a key concept in mathe- 
matical epidemiology or immunology. In fact, it de- 
termines whether the disease/virus dies out or persists. 
When R 0 > 1, the infection becomes chronic. In this 
case, it can be shown that the model supports a single 
"endemic" state. 
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